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ABSTRACT 

A large variety of jet properties can be described perturbatively by evolving 
the parton cascade down to small scales of order of a few 100 MeV. We dis- 
cuss two recent applications of this approach: 1. the soft limit of the particle 
spectrum which is nearly energy independent as expected for the soft gluon 
bremsstrahlung; further tests of this interpretation in e'^e~ and DIS are pre- 
sented; 2. the steep decrease of the distribution of rapidity gaps in e+e^ anni- 
hilation can be explained by the Sudakov formfactor for quark jets. 

1. Introduction 

An important goal in the study of high energy hard colhsion processes is the de- 
tailed understanding of the hadronic final state and its characteristic jet structure. 
There are two mayor fields of interest, one concerning the predictions in QCD per- 
turbation theory, mainly on the jet final states, the second concerning the colour 
confinement phenomena, how the partonic cascade evolves into the final state of 
hadrons. 

In the popular models for particle production in hard collision processes it is 
assumed that first the primarely produced partons evolve by bremsstrahlung processes 
according to perturbation theory into jets of partons until a characteristic resolution 
scale of Qo ~ 1 GeV is reached. Thereafter, non-perturbative processes take over and 
the final hadronic particles, often through intermediate resonances, are produced, for 
example, by a string mechanismS or through cluster formation.!! These models have 
been proven to be very successful in the description of the experimental data over 
the years. On the other hand, they are only accessible through a Monte Carlo code 
and involve quite a number of a priori unknown parameters in the description of the 
hadronization phase. 

We consider here another approach, based on the concept of "Local Parton Hadron 
Duality" (LPHD).i The parton cascade is evolved further down to a scale of about 
Qo ~ 250 MeV where Qo is now the cut-off in the transverse momentum. Then one 
compares directly an observed quantity with the corresponding parton level calcula- 
tion without any hadronization correction. 
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An immediate benefit of tliis approacli are tlie analytical results in many cases 
which provide a deeper understanding of how the theory is actually tested by a 
particular measurement. For example, the scaling properties and the systematics of 
their violations can only be formulated in a meaningful way using analytical formulae. 
Some results can be expressed in compact form by simple QCD numbers. 

Theoretical calculations can be carried out in the simplest case in the Dou- 
ble Log Approximation (DLA),Bi which provides the high energy limit, or in the 
Modified Leading Log Approximation (MLLA)i which includes finite energy correc- 
tions; they are usually essential to obtain quantitative agreement with experiment 
at present energies. The original success was in the description of the inclusive 
hadron distributions,! thereafter it has been applied to a large variety of different 
observables.0 In recent years more applications have been worked out and the ex- 
perimental results from LEP, HERA and TEVATRON gave further support to this 
approach.i 

The duality picture is clearly rather bold and it is not clear a priori for which ob- 
servables and in which kinematic ranges it applies. This has to be clarified ultimately 
by experiment. It is obvious that this model cannot compete with the standard 
hadronization models in the description of the various details of the final state like 
the production of different particle species or resonances. It has so far been applied 
successfully for suitably averaged inclusive quantities. 

The main virtue of this approach is its intrinsic simplicity which involves only two 
essential parameters, namely the QCD scale A and the non-perturbative mass scale 
Qo- III some cases the absolute normalization is taken as additional parameter; on the 
other hand, in a recent calculation of multiplicities beyond MLLA it was suggested 
that one parton corresponds to one hadron in this duality picture.! 

The hope is that with increased phenomenological knowledge about successes 
and limitations of this approach our theoretical understanding of the confinement 
mechanism will improve. In this presentation we dicuss two recent applications, one 
concerning the soft limit of the energy spectrum0 and the other one the rate for 
events with large rapidity gaps,0 i.e. events in a quasi exclusive limit. 

2. Main Ingredients of the Analytical Calculations 

In this section we give a short overview of the analytical treatment of multiparticle 
processes. This involves the matrix elements for the primary hard subprocess and the 
evolution of the primary partons into jets by bremsstrahlung cascades. The gluon 
bremsstrahlung at small angles d with energy E off the primary parton of type A 
[A = q, g) with jet momentum Pj is given by 

duA = J^^P^)^^, 7o(Px) = = ^^^^y >Qo 
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where E^, (3'^ = 4Nc/b, b = (llA^c- — 2nf)/3 with Nc, nj the numbers of 

colours and flavours, also Cg = Nq, Cq = 4/3. Inside the cascade the soft gluons are 
coherently produced from all harder partons. For azimuthally averaged quantities 
the consequences of the coherence effect can be taken into account by the angular 
ordering prescription^ which requires the angles of subsequent gluon emissions to be 
in decreasing order. 

The multiparticle properties of the jet can be discussed conveniently by using the 
generating functional Zj^{Pj,Q]u{p)). Here Pj and denote the initial parton mo- 
mentum and opening angle of the jet, and u{p) is a profile function for particle momen- 
tum p. The functional is constructed from all the exclusive final states. Then the in- 
clusive densities can be obtained by functional differentiation with respect to the pro- 
file function u{p), for example the one particle density by p^^\p) = 5Z{u} /5u{p) \u=i- 
The properties of these densities can be derived from the evolution equation for Z 
which relates the functional at scales Pj, G to the one at lower scales according to the 
"decay" A —>■ BC. In MLLA accuracy this evolution equation is given by0 
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dlnQ ' 2 ^ Jo 

x^^<l>Tiz)[Zs{zP„Q)Zc{{l-z)P„Q) - Za{P„Q)] (2) 

where denotes the DGLAP splitting functions. The initial condition for the 

parton evolution is given by 

ZA{Pj,&;{u})\p^e=Qo = UA{p = Pj), (3) 

i.e. at threshold there is only the primary parton. 

From (2) one can obtain the evolution equations for particle densities by appro- 
priate functional differentiation. Therefore this equation is the basic tool for deriving 
the multiparticle properties of a jet analytically. A more detailed discussion on the 



approximations can be found elsewhere. E 
3. The Soft Limit of the Particle Spectrum 



3.1. Theoretical Predictions Confronted with Experiment 

An essential ingredient of the LPHD approach is the evolution of the parton 
cascade towards small scales. The normalization of all quantities is given by the initial 
condition @. Therefore, one has to assume that the perturbative evolution can be 
applied at small cms energies as well. Alternatively, one could start the evolution 
from a larger cms energy whenever the perturbative QCD is considered trustworthy, 
but then one had to introduce unknown parameters for initialization instead of (^) 
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and a lot of the predictive power of this approach would disappear. Indeed, applying 
the MLLA calculations^ the moments of the energy spectra are well described down 
to such low energies as y/s = 3 GeV using the initial condition (^.0 

Of special interest is the behaviour of the soft end of the particle energy spec- 
trum. The coherence of the soft gluon emission from all harder partons forbids the 
multiplication of the soft particled^l and the spectrum in the variable ^ = \n{Pj/E) 
obtains a "hump-backed" shape.B0 

This problem has been studied recently in more detail.0 The parton density in 
the variable ^ = ln{Pj/E) in the DLA is given in lowest order of by 

duAjE^P,) _ Ca 2. HE /A) . . 

d\nE Nc^ ln(Qo/A) ^ ' 

This single bremsstrahlung contribution is independent of the jet energy Pj] the higher 
order contributions do depend on Pj but they are neghgable in the soft limit E Qq. 
Also the MLLA correction leaves this limit unaltered: the energy conservation effects 
and large z corrections from the splitting functions which make up the differences 
between DLA and MLLA can be neglected. If the perturbation theory and LPHD 
are valid towards such low particle energies E one expects then also that the hadron 
spectrum becomes independent of the cms energy in the soft limit. To exhibit the 
energy dependence it is convenient to consider the invariant density Edn/d^p which 
remains finite for small momenta." Then we expect an energy independent density Jq 
of hadrons in the limit where the particle momentum p, or alternatively rapidity y 
and transverse momentum k±, become small: 

/o= lim E^ = -limE^. (5) 

y^o,PT-*o d'^p 2 p^o d-^p 

The factor | takes into account that both hemispheres are included in the limit p 0. 

This prediction is a direct consequence of the coherence of the soft gluon emission: 
the emission rate for the gluon of large wavelength does not depend on the details 
of the jet evolution at smaller distances; it is essentially determined by the colour 
charge of the hard initial partons. 

In Fig. 1 we show the experimental results on the invariant density of charged 
particles for cms energies from 3 to 130 GeV in e~^e~ annihilation. An approximate 
energy independence of the low momentum particle density (within about 20%) is 
indeed observed; the same is true for identified particles vr, K and p.0 The curves 
in Fig. 1 represent the MLLA results. At very low energies where E k, Qq the 
QCD perturbative results have to be supplemented by a kinematical relation between 
parton and hadron spectra taking into account their different masses (for more details, 
see Ref. 10). The theoretical curves show the approach to the scaling limit and 

"Another alternative is the non- invariant spectrum dn/d'^pEi which has the advantage of being in- 
dependent of the particle mass. 
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Fig. 1. Invariant density EdnjdPp of charged particles in e^e annihilation as a function of the 
particle energy E = ^Jp^ + Qo *3o — 270 MeV. Experimental data at various cms energies from 
SLAC, TASSO and TOPAZ Collaborations, LEP-1 and LEP-1.5 are compared to MLLA predictions 
(taken from Ref. 10). 



describe well the different slopes at larger particle energies. An important role here 
is played by the running a,^ which provides the strong rise towards small energies for 
E <1 GeV, for fixed as this rise would be much weaker.lli0 

Recently inclusive spectra from HERA in the current fragmentation region in the 
Breit frame became available.lli The data are found to be rather similar in shape to 
the observation in e"'"e~annihilation at comparable scales. Especially, the approximate 
energy independence of the soft part of the spectrum is verified in the region 12 < 

< 100 GeV^. This supports the universality of the low energy phenomena. 



3.2. Further Tests of the Perturbative Picture 

One may suspect that the universal, approximately energy independent soft limit 
of the invariant spectrum could be some general hadronization phenomenon not re- 
lated to the coherence properties of perturbative QCD. It is therefore very important 
to further investigate the perturbative origin of the effect. Whereas there is no param- 
eter free prediction of the absolute size of the density Iq, the perturbative approach 
predicts its dependence on the colour charge of the primary parton as in (^). Espe- 
cially, we expect for the ratio of soft particle densities in a gg system to the one in a 
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qq system 




Such a ratio is expected asymptotically for the total multiplicities of particles,c3 but 
at present energies there are large corrections. On the other hand, such a large ra- 
tio may nevertheless appear in case of the soft particles where the DLA asymptotic 
behaviour is valid to a good approximation. The preparation of gg final states is not 
easy; but there are two lines of approach which can test the hypothesis (H), the 
first at a more qualitative, the second at a quantitative level. 

Study of the central rapidity plateau at small p± 

Let's first consider e"'"e~annihilation. The rate Iq in is defined in the cms. If we 
perform a Lorentz transformation in direction of one primary parton, we arrive at a 
frame with unequal momenta of q and q, say at rapidity y from the cms. As the soft 
radiation does not depend on the primary energy, the density Io{y) in the frame at y 
has not changed, so Io{y) should be constant away from the kinematical boundaries. 
Instead of the density at p± = we can also consider the density in an interval at 
small p±. This picture is a bit idealistic, in that the primary parton direction is not 
known and should be replaced by a suitable jet-axis, then Io{y) will not be exactly 
constant. 

Next we consider collisions with incoming hadrons or/and photons. The interac- 
tion may proceed through quark or gluon exchange where we have to assume some 
semihard momentum transfer of, say, at least 1 GeV — the soft processes may follow 
different rules. In case of a semihard quark or gluon exchange the outgoing partons 
will be either colour triplet or octet sources of radiation respectively, so in the soft 
limit we expect the ratio as in (^) for both cases. 

Let's now consider specifically the DIS process studied at HERA and construct 
Io{y) as above. Analogous results are obtained for 77 collisions. In the Breit frame we 
have a quark in the current direction recoiling against another colour triplet system 
in proton direction. So in this frame we have Io{yBreit) = Iq, the density as in 
e"'"e~annihilation. However, the same DIS process can also be initiated by other 
subprocesses, for example, the photon gluon fusion. This process can be viewed best 
in a frame closer to the proton: a qq pair moves in current direction and the remainder 
in proton direction. This process is mediated by gluon exchange, so the colour octet 
sources yield a soft particle rate Io{y) = |/o in such a frame. The overall picture 
then looks as in Fig. 2. 

As already mentioned the exchanged parton should be sufficiently hard. The 
gluon virtuality cannot be simply enforced by the external conditions but it should 
become larger with increasing Q^. We have no explicit calculation for this effect in 
the moment, so the conditions are formulated at a qualitative level. However, there 
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Fig. 2. Schematic view of particle density at small p± as a function of rapidity y in the cms: in 
frames at rapidity y with dominant gluon exchange the soft particle density Iq (y) is 9/4 times larger 
than in the frames with dominant quark exchange (the Breit frame, for example). 



are quantitative asymptotic expectations and predictions which could be falsified. We 
consider the distribution of particles in rapidity and at small p±, at best in the Breit 
frame. An increase of should make the small k± current plateau longer but not 
higher than Iq. At the same time the plateau in proton direction should be rising 
and approach for large and a limit, the "gluon plateau' 
larger than the "quark plateau" Jq. 



Iq, about two times 



Soft radiation perpendicular to the production plane in multi-jet events 
The collinear parton configurations just discussed correspond to quark or gluon ex- 
changes with the soft radiation expected in the ratio (|^). We consider now the con- 
figurations not exactly collinear but with moderate transverse momentum, so that 
the corresponding parton jets can be resolved. In that case we consider the soft 
particle production in direction perpendicular to the production plane which can be 
analytically calculated in the soft region taking into account the cut-off Qq. 

The simplest case is the process e'^e^ qqg. The soft radiation in this process 
has been analysed already some time agoi and applied to the case of multiplicity 
flow in the production plane. The spectrum perpendicular to the plane (say, in a 
small cone) is conveniently normalized to the corresponding flow in a 2 jet qq event 
by 



RUp) = 



dNl/dQridp 
dNf/dVLfidp 



(7) 



for a general process a. For the case of qqg events one flnds for the soft bremsstrahlung 
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in order 



RTiv) 



Nc 



2 - cos Qgg - COS Ogg - -^(1 - COS Qgg) 



where Qij are the angles between the jets in the cms. This formula yields the two 
collinear limits = 1 for collinear or soft gluons and R_[_ = | for antiparallel 
quarks recoiling against the gluon. In this approximation the shape of the momentum 
spectrum does not depend on the inter-jet angles but the absolute magnitude does. 

The analogous analysis can be carried out for photo-production of dijets in jp or 
77 collisions and also in pp collisions.lii One can distinguish between direct and re- 
solved pro duct ion. ii'ii The former involves the soft radiation from a qqg system, very 
much like in e~^e~ qqg, and for the dominant small angle scattering corresponds 
to quark exchange, whereas the latter at small angles corresponds to gluon exchange. 
One expects again the ratio (||) of the soft particle production in both processes at 
small scattering angles. 

It will be interesting to find out whether the particle density varies with the angles 
as predicted by (H) and occurs in the predicted ratios in different processes even down 
to low momenta of a few hundred MeV. This would provide a strong argument in 
favour of the perturbative QCD interpretation of the soft particle production. 



4. Large Rapidity Gaps in e+e Annihilation 

In the previous example of inclusive particle spectra the partons and hadrons are 
typically close in phase space and the hadronization requires only small rearrange- 
ments of momenta. If during the evolution of the parton cascade in e"'"e~ annihilation 
a large rapidity gap occurs, one may expect the colour confinement forces at large dis- 
tances to fill the gap by hadrons in the hadronization phase. Then there is no close 
relation between parton and hadron final states. The study of large rapidity gap 
events can give us therefore further insight into the colour neutralization mechanism. 

The interest in such events has been awakened recently by the findings of large 
rapidity gaps in dijet events at TEVATRONil and HERAii which are interpreted by 
the exchange of colour neutral objects. 

In e^e~annihilation — if one wants to create a large gap without colour exchange 
such as to avoid final state rearrangements by confinement forces — one is lead to hard 
processes of the type e~^e~ — >■ qqqq or qOQQ where the hadrons are produced from 
colour singlet low mass qq or gg clusters.B3'E3 As these processes in the perturbative 
analysis involve a highly virtual intermediate quark or gluon their rate is rather small 
though. One finds that the rate keeps decreasing with increasing gap size, contrary 
to the case of pp or ep collisions. 

Indeed, the recent measurement of rapidity gaps by the SLD coUaborationil has 
shown an unlimited decrease of the gap rate over five orders of magnitude (such a 
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decrease has also been seen at lower energies some time ag However, thegap 
rate in absolute terms exceeds the expectations from the above calculations0'0^by 
about two orders of magnitude and therefore, this type of hard colour neutralization 
cannot be the dominant process for the formation of large rapidity gap events. 

An alternative approach,0 in line with the previously discussed LPHD, considers 
a rather soft confinement mechanism, i.e. the hadron formation does not change the 
parton final state considerably, not even if a large rapidity gap is formed. Clearly, 
such an application of LPHD cannot be justified a priori. Rather one tries and looks 
how different the rates of large gaps really are in the parton and hadron cascades. 

One has to calculate the probability that no parton is emitted into a certain 
angular interval, say between 9i and Q2 (©i < ©2) where 9 is measured with respect 
to the initial parton direction. The rapidity for a massless parton is then obtained 
from y = — Intg^. This probability is given by the so-called Sudakov formfactor. 
For the actual problem of the angular ordered cascade with transverse momentum 
cut-off it has been calculated in connection with the multi-jet rates in the Durham/p^ 
algorithm .Bli 

The probability that no parton be emitted below the angle in the cascade with 
cut-off p± > Qo initiated by the primary parton A of momentum Pjis given by 

A^(P„e,go) = exp{-WAiPj,Q,Qo)) (9) 
wa{Pj,&,Qo) = [doo'f de'pA{uj',e'), (10) 

where Pa{^', ©') = duA/du'dQ' is the emission density of a gluon at energy uj' and 
angle 9'. The probability for a gap between the two angles 9i < 62 < f is then 
given by the ratio of two Sudakov formfactors 

/A(ei, 82) = AA(e2)/AA(ei) = exp (-^^(62) + WAiQi)) . (11) 

In the simplest approximation, the DLA, one finds 

wa{Y, A) = ^f^'HY + A) In ^ - Y} (12) 

In this approximation the jets evolve independently in the two hemispheres and 
the respective gap probabilities factorize. In particular, for the symmetric gap in 
e"'"e~ annihilation one obtains accordingly 

/^(Gg) = e-2(--4(i)-"'^(®«)). (13) 

A result with improved accuracy is obtained by taking into account the exact form of 
the matrix element for e~^e~ qqg with numerical evaluation of the integral (|1OD.0 
An analytical approximation in MLLA can be obtained by simplifying the integral 
over the nonsingular parts of the splitting functions.lll0 
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Fig. 3. Fraction of rapidity gap events in e+e" annihilation as a function of the fuh width of the 
symmetric gap. The data points refer to the measurement of gaps between charged particles (r- 
lepton events included) by the SLD Collaboration.c3 The dashed histogram shows the expectation 
from the JETSET El Monte Carlo without r-leptons. Also shown are the DLA predictions for the 
gaps of a parton cascade for two values of A = ln(Qo/A at Qq = 0.244 GeV. 



As the result of these calculations one finds the probability for the gap in the 
symmetric rapidity interval Ay ^ — 21n^. It shows an almost exponential decrease 
very much like the experimental data. The slope is not very different for DLA and 
MLLA calculations but strongly depends on the parameter A = ln((5o/A). 

The SLD dataS are shown in Fig. 3. They refer to rapidity gaps between charged 
particles and still include r-lepton events. The analytical calculations^ from the DLA 
are shown for two parameters of A at A = 0.244. An upper limit A < 0.1 has been 
derived from a moment analysis of the spectra.0 These calculations should rather be 
compared to data on gaps empty of any particles, not only of charged ones. Their 
distribution is expected to have a steeper slope than shown by the data in the figure. 
As an example the prediction for A = 0.05 is shown. A recent analysis of particle and 
jet multiplicities beyond has given the small value A = 0.015 ± 0.005. 

One can conclude that the distribution of gaps in a hadron cascade is well rep- 
resented by the one in the parton cascade, provided the cut-off Qq is taken small 
enough, i.e. very close to the QCD scale A in the present approximation scheme. If 
one took the larger cut-off Qo ~ 1 GeV the gap probability at Ay = 6 would be 
larger by 4 orders of magnitude! This suggests that the hadronization phase in a 
conventional hadronization model with shorter parton cascade has a similar effect as 
the evolution of the parton cascade from Qq = 1 GeV down to Qo = 0.25 GeV. This 
is not too surprising as it is in this last phase where the running increases most 
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Fig. 4. In the standard hadronization model (left) the parton cascade is followed by a hadronization 
phase. In the dual approach the hadronization phase is assumed to be represented by the parton 
cascade in an average sense after evolution towards a small cut-off scale Qo 

strongly and makes the suppression of radiation increasingly difficult. 

An interesting prediction from this approach concerns the distribution of rapidity 
gaps in single quark and gluon jets. As the exponent of the Sudakov formfactor is 
derived from the one-particle-density it is proportional again to the colour factor Ca, 
i.e. the slope of the gap distribution is correspondingly steeper (roughly twice) for 
gluon jets than for quark jets. This prediction could be checked in high p± jets at the 
TEVATRON or at HERA. 

5. Conclusions 

The description of detailed properties of the hadronic jets in hard processes in 
terms of a parton cascade assuming a duality between parton and hadron production 
continues to be successful. It is surprising that this duality works even for observables 
where a priori there is no good reason to believe in a pcrturbative approach. 

It appears that the evolution of the parton cascade at very low scales below 1 GeV 
resembles in its main inclusive aspects the hadronization phase of the phenomeno- 
logical models with all its resonance decays, so one can speak of a dual description 
of the hadronic phenomena after appropriate averaging (see Fig. 4). As the cut- 
off Qq is close to the QCD scale A within less than 10% the coupling ag becomes 
large in the last phase of the perturbative cascade. So it may not be unreasonable 
that the strong coupling perturbative regime corresponds to the hadronic phase with 
resonance production. 

The QCD expectation of the energy independence of the low momentum particle 
production from the soft gluon coherence is well met by the hadronic data in the 
full energy region explored, i.e. in e"'"e~annihilation from ADONE at 1.6 GeV up to 
LEP-1.5 and in DIS from = 12 GeV^ to = 100 GeV^. It wiU be interesting to 
test the predictions on the higher production rates of soft particles for other colour 
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configurations than qq which can be obtained in various processes. 

The parton cascade with low cut-off Qo and without hadronization also repro- 
duces the approximately exponential decrease of the rapidity gap distribution in 
e"'"e~annihilation. This suggests in particular that the confinement of colour in the 
jet evolution is a soft process. 

The correspondence between the parton and hadron final state according to these 
results is not like the conventional duality between a colour singlet qq and a hadronic 
cluster. Rather it is a correspondence between a coloured parton (typically a gluon) 
and a colourless hadron. These partons are taken at a particular scale Pj_ > Qo ^ A, 
so that partons below that scale are not resolved. These soft non-pcrturbativc sca- 
partons could take the important role in the confinement process without disturbing 
the momenta of the resolved perturbative partons. 

Meanwhile it seems to be a challenge to test further this dual picture of hadron 
formation at the limit of what is accessible perturbatively. 
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